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Abstract 

Non-spherical dynamical approximations and models for the gravitational collapse are used to extend 
the well-known Press &; Schechter (PS) approach, in order to determine analytical expressions for the mass 
function of cosmic structures. The problem is rigorously set up by considering the intrinsic Lagrangian nature 
of the mass function. The Lagrangian equations of motion of a cold and irrotational fluid in single-stream 
regime show that the shear, which is non-locally determined by all the matter field, is the quantity which 
characterizes non-spherical perturbations. The Zel'dovich approximation, being a self-consistent first-order 
Lagrangian and local one, is used as a suitable guide to develop realistic estimates of the collapse time of 
a mass clump, starting from the local initial values of density and shear. Both Zel'dovich-based ansatz and 
models and the homogeneous ellipsoidal model predict that more large-mass objects are expected to form 
than the usual PS relation. In particular, the homogeneous ellipsoid model is consistent at large masses 
with a Press &; Schechter mass function with a lower value of the 6c parameter, in the range 1.4-M.6. This 
gives a dynamical explanation of why lower 6c values have been found to fit the results of several N-body 
simulations. When more small-scale structure is present, highly non-linear dynamical effects can effectively 
slow down the collapse rate of a perturbation, increasing the effective value of 6c- This may have interesting 
consequences on the abundance of large-mass high-redshift objects. 

Subject headings: cosmology: theory — galaxies: clustering — large-scale structure of the universe 
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1 Introduction 



A fundamental feature of our observable Universe is the presence of collapsed clumps of matter, such as galaxies, 
groups and clusters of galaxies, whose mass distribution (mass function, hereafter MF), which has recently been 
observationally estimated (see, e.g. Ashman, Salucci &: Persic 1993, for galaxies; Pisani et al. 1992, for groups; 
Henry &: Arnaud 1991, for X-ray clusters; Bahcall &: Cen 1993, and Biviano et al. 1993, for optical clusters), 
is a relevant quantity for comparison with the predictions of cosmological models of structure formation. Much 
effort has been focused on the problem of finding the theoretical MF of a given model; this is not an easy task, 
as the equations of motion of matter in the collapsing phases are non-local and highly non-linear. 

A first attempt to estimate the MF of collapsed structures has been made by Press &: Schechter (1974, 
hereafter PS) by means of very simplified assumptions: given a smoothed initial gaussian density field in a flat 
Universe, every mass element is supposed to follow spherical collapse (e.g., Peebles 1980), according to which 
infinite density is reached when the density contrast 6, linearly extrapolated to the collapse time, is 1.69; in 
this context, the dynamical fate of a perturbation is dictated simply by its initial density contrast. It is then 
supposed that the collapsed element forms a structure with mass equal or greater than the mass associated to 
the smoothing volume adopted, so that merging of collapsing mass elements into larger ones is allowed. This 
approach has three main weak points: (i) the dynamical model of spherical collapse is unrealistic in general cases, 
(ii) The PS mass function fails to predict the collapse of clumps which are not overdense enough to collapse 
by themselves, but nonetheless are going to be accreted onto larger structures [cloud-in- cloud problem); then, 
to reach the correct normalization (all matter collapsed at some scale), one has to multiply the MF by a 
"fudge factor" 2. (iii) Every mass element evolves independently from the others, so that coagulation and 
fragmentation^ of collapsed clumps are neglected. 

Weak point (i) has not generally been faced: most of the results on the MF and on the abundance of 
structures, such as rich clusters of galaxies or quasars, rely on spherical collapse, which is asymptotically 
accurate only for structures coming from the rare highest peaks of the primordial density field (e.g. Bernardeau 
1994). Sometimes the homogeneous ellipsoid model has been used to describe supercluster dynamics (e.g. 
White &: Silk 1979, Hoffman 1986). More recently. Bond &: Myers (1993a, b) have used the ellipsoidal model to 
describe the collapse of peaks, while Bartlemann, Ehlers &: Shneider (1993) have used both ellipsoidal collapse 
and Zel'dovich approximation to estimate the collapse time of a mass element. 

On the other hand, points (ii) and (iii) have been faced by many authors. The normalization "fudge factor" 
2 has been justified by Peacock &: Heavens (1990) and Bond et al. (1991) by means of the excursion set theory; 
briefly, collapse prediction and mass assignment have to be done on the basis of the largest collapsing structure 
which includes the given initial point. However, the resulting MF is the same as the PS one only in the case 
of sharp ^space smoothing. Furthermore, coagulation and fragmentation of clumps have been described by 
means of Smoluchowsky equation (Silk &: White 1978; Lucchin 1988; Cavaliere, Colafrancesco &: Menci 1991; 
Cavaliere &: Menci 1994), with the result that at small scales, or in some special kind of environment (e.g. cD 
galaxies in rich clusters), the MF can lose memory of the cosmologically relevant initial conditions. Finally, 
Cavaliere &: Menci (1994), by means of Cayley trees formalism, have succeeded in unifying the excursion set 
and Smoluchowsky equation approaches. 

Full N-body calculations have surprisingly shown that the PS mass function gives reasonable fits, despite 
of its weak points (Efstathiou et al. 1988; Efstathiou &: Rees 1988; Carlberg &: Couchman 1989; Bond et al. 
1991; Bond &: Myers 1993b; Lacey &: Cole 1994; Klypin et al. 1994); however, the comparisons between PS and 
N-body MFs have been made only in limited ranges of masses, and some tuning of the MF parameters has been 
needed to achieve good fits. At variance with the previous authors, Jain &: Bertschinger (1993) have claimed 
that, for a standard cold dark matter spectrum, the PS mass function significantly underestimates the number 
of high-redshift objects, while it gives a reasonable fit at later times, consistently with their calculations on the 
non-linear growth of perturbations. They have ascribed the difference between their result and the previous 
ones to the improved resolution of their N-body P^M code. 

^Here two mass elements merge if they are part of a larger-scale collapsing clump, coagulate if their interaction is rather a 
two-body collision between already collapsed clumps; coagulation and fragmentation are responsible for that part of the merging 
histories which can not be predicted from local collapse and mass assignment criteria. 
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The PS approach, originally limited to an Einstein-de Sitter Universe, has been extended to the more 
general cases of open Universes with non-vanishing cosmological constant (Lilje 1992, and references therein), 
and of non-Gaussian statistics (Lucchin &: Matarrese 1988). Other approaches have been used to determine an 
analytical expression for the MF. Within the peak paradigm, according to which structures form in the peaks 
of the initial density field, the density of peaks (Bardeen et al. 1986) can be used to construct a MF; however, 
it is not clear which mass has to be assigned to a peak and, moreover, the peak paradigm has been put under 
discussion by recent N-body works (Katz, Quinn &: Gelb 1993; van de Weygaert &: Babul 1994; but see, for a 
different view. Bond &: Myers 1993a, b). Finally, Cavaliere, Colafrancesco &: Scaramella (1991) have developed 
a dynamical approach to the mass function, explicitly introducing time-scales for creation and destruction of 
clumps. 

To deepen the knowledge on the MF it is necessary to consider more realistic models for the collapse of a mass 
element than the spherical one. In this paper I consider dynamical models and approximations for the non-linear 
evolution of perturbations which are able to predict the fate of a mass element from local initial conditions; in 
this case it is possible to extende the original PS approach, which seems to offer the best compromise between 
simplicity and success. In §2 the intrinsic Lagrangian nature of the MF is stressed, so that the Zel'dovich 
approximation, which is local, provides a natural and useful guide to estimate the collapse time of a mass 
element. In §3 a correction to the PS mass function is defined and some ansdtz for the collapse time, based on 
the Zel'dovich approximation, are introduced and their consequences on the MF are examined; then, two more 
realistic models are analyzed: one based on the Zel'dovich approximation, and the other on the homogeneous 
ellipsoid model. In §4 the physical differences between the usual PS mass function and the new ones here 
proposed are outlined and discussed, and in §5 the results are summarized and final conclusions are drawn. 



2 The Mass Function As A Lagrangian Quantity 

To assess in a rigorous way the problem of finding an approximation for the collapse time, it is useful to start 
with the equations for the evolution of a cosmological fluid. What is relevant to our purpose is the density of a 
given mass element — the Lagrangian density — and, in particular, the instant at which it diverges, which is 
the collapse time. Then, the MF can be regarded as an intrinsically Lagrangian quantity. It is useful to restrict 
ourselves to a matter-dominated flat Universe with vanishing cosmological constant, filled with pressureless and 
non-vortical matter; moreover, the initial field is assumed to be linear, Gaussian, and growing-mode dominated, 
and the dynamical description is limited to the single-stream regime (no orbit crossing having occurred yet). 
These assumptions, which are shared by the usual PS approach, can be relaxed to more general cosmologies, 
retaining the validity of the qualitative features of the present results. 

A convenient form for the Lagrangian Newtonian evolution equations of a cosmological fluid can be obtained 
(Gurbatov, Saichev &: Shandarin 1989; Matarrese et al. 1992) by using the scale factor a (which is also the 
growth factor of the linear growing mode) as a time variable (a is normalized to be unity at the present time). 
In this way, being x the Eulerian coordinate, the velocity of the fluid element is u = dx/da = x = w/a{da/dt) 
(v is the usual peculiar velocity), d/da = d/da + u- Vx is the Lagrangian time derivative, the density contrast 
6{q) is defined as usual as [£'(q) — and the potential (p is defined as ^ = 3tQ/2aQ(f), where (f) is the usual 

gravitational potential. With these definitions, the Euler, continuity and Poisson equations become: 

3 3 

u+— u = - — Vx¥' (1) 
2a 2a 

5+(l + 5)Vx-u = (2) 

Vi<p = - (3) 
a 

These equations are valid as long as the fluid is in single-stream regime. 

The Eulerian gradients of the components of the peculiar velocity u can be decomposed, without loss of 
generality, into the expansion scalar i?, the shear tensor S and the vorticity tensor w: 
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ViUj = -dSij + + uiij (4) 

(see, e.g., Ellis, 1971). i? is simply the divergence of u, while S,j and w,j are the traceless symmetric and 
antisymmetric parts of V, Wj . w,j is of course null for an irrotational fluid. Whit these kinematical quantities it 
is possible to find the evolution equation for the density contrast: 

where = and = u>ijU>ij. In the case S = w = 0, Eq. (5) describes the evolution of a spherical 

perturbation, and its solutions are the well-known ones of the spherical collapse model used in the PS mass 
function. It is clear that, while the w = hypothesis is quite reasonable in a cosmological context, at least before 
orbit crossing (e.g. Peebles 1980), the hypothesis of vanishing shear is unrealistic for a general perturbation. It 
is to be noted that the shear-dependent term in Eq (5), S^, is positive-definite, so that, as long as the fluid is 
irrotational, the growth rate of the density contrast of a mass element is always enhanced by the shear. As a 
consequence, a general mass element always collapses faster than an equally overdense spherical one. This fact, 
already noted by Hoffman (1986), has been put in form of theorem by Bertschinger &: Jain (1994). Moreover, 
the shear term is the one which introduces the non-locality of Poisson equation into Eq. (5); indeed, the 
evolution equation for S,j is : 

±ij + + + ^S.j - ^^^6ij = (viVj<p - ^V\6ij^ , (6) 

where the source term between parentheses on the rhs is the tidal field tensor, which is non-local (e.g., Kofman 
&: Pogosyan 1994). Non-local estimates of the collapse time make the PS approach inappropriate and introduce 
strong computational difficulties; then, in order to study the effect of non-spherical collapse with a PS-like 
approach, only local estimates of the collapse time will be considered in the present work. 

A number of recent papers have focused on the fact that the GR Lagrangian evolution equations of a 
pressureless irrotational fluid are local if the so-called magnetic part of the Weil tensor is neglected (Barnes 
&: Rowlingson 1989; Matarrese, Pantano &: Saez 1993; Bertschinger &: Jain 1994). As the Newtonian limit of 
that tensor is null, it can seem natural to extrapolate this approximation to that limit; however, non-null post- 
Newtonian terms from enter in the Newtonian equations of motion introducing non-locality (Matarrese, 
Pantano &: Saez 1994; Kofman &: Pogosyan 1994), so that this = approximation is not valid in general 
but only for some special symmetries (Bertschinger &: Jain 1994). 

Anyway, it is possible to find well-founded approximations which give a local description of the evolution 
of a mass clump. In particular, the most famous approximation which can be introduced in the Lagrangian 
framework, the Zel'dovich one (Zel'dovich 1970, hereafter ZEL; see also Shandarin &: Zel'dovich 1989) is local. 
Let q and x be the Lagrangian and Eulerian coordinates of a fluid element, with x(q, ao) = q; the trajectory 
of the element is: 



x(q,a) = q+*(q,a), (7) 

where ^ is called displacement field. ZEL consists in expressing the latter as the product of a universal time 
factor b[a) and a constant field u(q): 



x(q,a) = q + 6(a)u(q); (8) 

here b is the growth factor of the linear growing mode, b[a) = a — ao (in a flat Universe), and the field u(q) is the 
initial peculiar velocity of the fluid element. We will refer to Eq. (8) as ZEL mapping. According to it, the fluid 
elements maintain their initial peculiar velocity (li = 0), performing straight paths in the Eulerian comoving 
space. In this way, all physical quantities can be derived from the local values of the initial deformation tensor, 
i.e. the Lagrangian derivatives of the peculiar velocity field, dui/dqj. 
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At a certain time fluid elements coming from different Lagrangian positions get to the same Eulerian point 
(orbit crossing or shell crossing); there caustics (i.e. regions of inflnite density) form, and ZEL mapping becomes 
multi- valued, i.e. multi-stream regions form. In this case ZEL predicts the particles just to continue their straight 
walks, while in the real collapse the particles will form a structure. I have chosen to identify the collapse time 
with the flrst shell-crossing, as it is the instant at which the density of the mass element diverges; moreover, 
shock waves in the baryonic component are predicted to form at that moment (Zel'dovich 1970). This point is 
discussed in a deeper way at the end of §4. 

ZEL can be also obtained as the flrst-order term of the expansion of the displacement fleld ^ in terms 
proportional to 6(a)" (e.g., Buchert 1994 and references therein): it is a true self-consistent flrst-order Lagrangian 
approximation. The fact that the series expansion of the displacement fleld and not of the density contrast 
is made implies that ZEL, together with higher-order Lagrangian approximations, is not a priori limited to 
5 <C 1; as a matter of fact, ZEL is an exact solution in ID, is accurate up to 5 ~ 1, and gives answers in 
the right direction also in the highly non-linear regime up to shell crossing (e.g., it does predict caustics to 
form, while the usual Eulerian linear approximation does not). Higher-order terms of the displacement fleld 
can be calculated, and their convergence toward the exact solution is assured up to orbit crossing. However, 
higher-order Lagrangian approximations are non-local, so they will not be considered here. 

In the realistic case of matter flelds which have power on all scales, so that multi-stream regions can always 
be found at small scales, to apply ZEL it is flrst necessary to smooth the initial density fleld, or, in other words, 
to truncate the power spectrum at some wavenumber (Coles, Melott and Shandarin 1993). This truncated 
version of ZEL has been found, by means of cross-correlation tests, to predict in a satisfactory way the relevant 
large-scale features of N-body simulations for spectral indexes n not greater than 1 (Coles et al. 1993; Melott, 
Pellman &: Shandarin 1993; see also Melott 1994 and references therein). Moreover, the less small-scale structure 
is present, the best truncated ZEL works; gaussian smoothing optimizes its performances. Furthermore, peaks 
in ZEL-evolved flelds have been shown to describe fairly well cluster correlations (Mann, Heavens &: Peacock 
1993; Borgani, Coles &: Moscardini 1993). Other authors (Grinstein &: Wise 1993; Bernardeau et al. 1993, 
Munshi &: Starobinsky 1993, Munshi, Sahni &: Starobinsky 1994), by means of perturbative calculations, have 
shown that ZEL underestimates the higher moments of the evolved density fleld (i.e. the skewness, kurtosis 
etc.), but nonetheless it works better than other approximations schemes. 

Other approximations for the non-linear evolution of density perturbations have been proposed, but they 
are generally non-local. For the adhesion model (see, e.g., Shandarin &: Zel'dovich 1989), in which particles 
whose trajectories try to intersect simply stick together, Vergassola et al. (1994) have analytically calculated the 
MF with the aid of several mathematical tools borrowed from diffusion theory, flnding a large-mass behaviour 
similar to the PS one, but a different small-mass slope. Other approximation schemes, such as frozen flow 
(Matarrese et al. 1992) and frozen potential (Brainerd, Sherrer &: Villumsen 1993; Bagla &: Padmanabhan 
1994) approximations, are intrinsically Eulerian, so that particle trajectories have to be found to get any 
dynamical prediction for a given fluid element. In the following, the self-consistent ZEL approximation will be 
considered as a guide to realistic estimates of the collapse time. 

3 Local Approximations for the Collapse Time 

3.1 A PS-like ZEL-based Approach to the MF 

Given a local estimate of the collapse time and the statistical probability distribution function (PDF) of the 
initial conditions, the PS approach gives a recipe to construct an analytical expression for the MF. As a flrst 
step, ZEL is used to approximate the dynamics of a mass element. It has been shown before that the Lagrangian 
peculiar velocity of a fluid element does not change in time, so that: 

u(q) = -Vq^o, (9) 

where ipo is the initial peculiar gravitational potential. If we call Ai, A2 and A3 the three eigenvalues of the 
deformation tensor dui/dqj (or, equivalently, of the tensor —d^ipo/dqidqj), assuming Ai > A2 > A3, the Jacobian 
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matrix Jij of the ZEL mapping, Eq. (8), can be written as: 



dx' du' 

Jij = -^ = 6ij + h{a)—^ = diag (1 + bXi, 1 + 6A2, 1 + 6A3) . (10) 
aqj aqj 

Equation (9) (9) and Poisson Eq. (3) imply Ai + A2 + A3 = —6o/a,o', then the three eigenvalues A,- give a full 
description of the initial conditions for the local evolution of a density perturbation. However, it is worth noting 
that these A,- eigenvalues are not local, in the sense that they are non-locally determined by the initial tidal field. 
This means that, at variance with the spherical top-hat model, which assumes that the Universe is homogeneous 
outside the perturbation, our mass elements are extracted from a perturbed density field and their dynamics is 
influenced, at least through the initial conditions, by the rest of the Universe. 

It is very easy at this point to explicitly calculate the expansion and the shear of a mass element: 

1? = (/Xl + /X2 + Ms) 

S.J = diag (/xi - l/3i?, fi2 - l/3i?, fi3 - l/3i?) ^ ^ 

where /x,- = A,/(l + aA,), i=l,2,3 (note that it has been assumed b[a) = a — ao ~ a). 

There are two different ways to deflne a density with ZEL (Shandarin, Doroshkevich &: Zel'dovich 1983). 
First, the continuity Eq. (2) can be used to construct a density which is consistent with mass conservation; 
recalling that the Jacobian determinant represents the volume of a mass element, it is easy to obtain: 

This continuity density 6c is not consistent with momentum conservation; we can deflne a dynamical density 6d 
by means of the Poisson Eq. (3), and by means of Eq. (9): 6d = — ai?. In this way: 

6d = -a{fjLi + fjL2 + fJ-a)- (13) 

6d is consistent with momentum conservation but not with mass conservation. All these quantities, i.e. the 
expansion, the shear and the two densities deflned above, diverge when 1 + aAa = 0, where A3 is the smallest of 
the A,-. Then the collapse time ac can be easily written as: 

ac = (14) 

Given this local estimate of the collapse time, to calculate the MF we need to know the PDF of the initial 
conditions, i.e. of the three eigenvalues of the deformation tensor A,-. This has been calculated for a Gaussian 
fleld by Doroshkevich (1970): 

675a/5 / 3 15 \ 

V{Xi,X2, A3) = — ^ exp -si + —S2 (Ai - A2)(A2 - A3)(Ai - A3); (15) 

oTTcr y (7 Zc J 

where cr is the usual mass variance linearly extrapolated to the present time, si = Ai + A2 + A3 and S2 = 
A1A2 + A2A3 + A1A3. It is convenient to perform the following change of variables: 

= — Ai — A2 — A3 

= A1-A2 (16) 
= A2 — A3 

where 6 is the density contrast linearly extrapolated to the present, and varies from -co to co. The variables 
X and y have been deflned in order that their ranges of variability are decoupled: the two inequalities A2 < A3 
and Ai < A2 implicate the conditions a; > and y > 0. They are simply related to the shear tensor: it is easy 
to show that the initial eigenvalues of S,j are 2x/3 + y/3, —x/3 + y/3 and —x/3 — 2y/3. Then, if a; = and 
y = the shear is null and the collapse is spherical. 
With the transformation (16) the PDF becomes: 
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V{6,x,y) = ^^exp I --^ ) x ^ W exp I --^(x^ + + y^) ) xy{x + y)=P{6) xP{x,y) (17) 



Equation (17) is just the product of two normalized distribution: the one for 5, 'P(5), which is the usual Gaussian 
distribution, and the one for x and y, V{x,y). 

With these new variables the ZEL collapse time becomes: 

ac = 7 . (18) 

6 + x + 2y ^ ' 

It is useful to define a function 6c{x,y) as the density contrast needed by a perturbation to collapse at ac = 1: 



6c : ac{6c{x,y),x,y) = 1. 



(19) 



In the ZEL case this is: 



6c{x,y) = 3-{x + 2y). (20) 

This 6c{x,y) function corresponds to the PS parameter Sc'- much more dynamical information is present now. 
As expected, the shear lowers the value of 5c, i.e. it helps the collapse: with an adequate shear even underdense 
perturbations can collapse. Moreover, rare density peaks, which are great overdensities with respect to the 
density variance, are in general characterized by a low shear, so that their collapse is nearly spherical; this is in 
agreement with Bernardeau (1994). 

In the case of spherical collapse, ZEL predicts 6c = 3, at variance with the exact value 1.69. This is not 
unexpected, as ZEL is known to underestimate the true collapse, especially in the spherical case. For the 
moment, let's assume that an "exact" 6c function is given, in the sense that it correctly gives, as a function of 
the local variables x and y, the initial density contrast necessary to make a perturbation collapse now. It is 
possible to write it as: 



^c{x,y) = 6o- f{x,y) 



(21) 



where 6o is the 1.69 spherical value, and f{x,y) is a positive definite function, as the shear accelerates the 
collapse, with /(O, 0) = 0. 

As in the PS approach, it is assumed that a collapsing mass element, whose associated mass is a function of 
the mass variance a through the smoothing scale R, M = M[R[(t)), becomes part of a structure whose mass 
is greater than M . The cumulative MF, i.e. the fraction of matter contained in structures with mass greater 
than M, is then simply the integral of the PDF, V{6,x,y), over all the initial conditions which predict the 
perturbation to have collapsed within the present time, i.e.: 



/lOO /lOO /lOO 

F{>M)= dx dy d6V{6,x,y) 

Jo Jo J6c{x,y) 



(22) 



The mass function is then: 



N{M)dM 



d ^, ,,,dM 
-Q F{> M) 



dF da- dM 



(23) 



^ da dM M 

In this equation, the term which contains the dynamical information is dF/da. The integral of the PDF 
over 6, being V{6) Gaussian, gives simply a term proportional to erfc(5c(a;, i/)/\/2o'), where erfc is the usual 
complementary error function. Integrating Eq. (22) in the variables X = x/a and Y = y/a, the only a 
dependence of the cumulative MF F is in the erfc term. We have: 



da \ V2cr 



27r(7' 



■ exp 



20-2 



aX 



d{aX) 



aY 



d{aY) 



(24) 
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In the PS case the term dF/da is: 



PS sPUct"^ V 2(7 
With Eqs. (17, 21, 22, 24, 25), the dynamical term dF/da can be written as: 



'° expf-i|). (25) 



dF _ dF 
da da 



X I{a), (26) 



PS 



where X((7) is a correction factor to the PS mass function which gives the contribution to a given mass scale of 
non-spherically collapsing clumps: 



/»00 /»00 / -t C \ 

X(<7) = y^ dXj^ dYV{X,Y)expl- — f{aX,aY) + -^f{aX,aY)\ 



^g{aX,aY) 

00 



(27) 



where g[x, y) = f — xdf/dx — ydf/dy. This expression is very general: given any local estimate of the collapse 
time, a function 6c can be constructed by means of Eq. (19); then the corresponding correction factor I{a) to 
the PS mass function is given by Eq. (27). 

A remark on the normalization: the fact that the PS normalization is wrong by a factor 2=1/0.5 is simply 
due to the fact that only overdense regions, which account for 0.5 of the mass, are predicted by the spherical 
model to collapse at some time. In the case of ZEL approximation, all the mass elements with negative As-values 
are going to collapse at some time, and this happens for 92% of the mass (Doroshkevich 1970). It could then be 
attempted to get the right normalization by multiplying the ZEL-based MF by a factor 1/0.92~1.09. However, 
if excursion set theory has justified the factor 1/0.5 for the PS mass function, it is not guaranteed that it would 
justify an analogous one 1/0.92 for the ZEL case. As a matter of fact, we are interested in the large-mass part 
of the MF, where a PS-like approach is more believable, and whose structures collapse in a quasi-spherical way, 
the case in which the fudge factor 2 is justified. Then the MF has been normalized so as to recover the correctly 
normalized PS one in the spherical case, which is 6c = 6o, = 1, i-e. the PS and corrected MFs have been 

multiplied both by the same fudge factor 2. This causes the corrected MF to have a wrong normalization by a 
factor 1.84, which is due presumably to an overcounting of objects which collapse in a definitely non-spherical 
way, i.e. small-mass objects. 

3.2 Realistic Estimates of the Collapse Time 

As stated before, the plain ZEL prediction for the 6c function does not reduce to the correct spherical value in the 
zero-shear case. Fig. 1 shows the ZEL 6c curve, Eq. (20) (dashed line), as a function of the variable z = x + 2y. 
A simple way to force 6c to reach 5o=l-69 at z=0 is to join it with the (0,5o) point by means of a horizontal 
branch: 6c = min{5o, 3 — z}. Physically, the mass element is supposed to follow spherical collapse, unless plain 
ZEL prediction allows it to collapse faster. Other possible 6c curves, which have the same normalization, can 
be obtained by shifting the plain ZEL prediction horizontally, always connecting it to 6o with a horizontal line 
(dotted lines in Fig. 1). This can be expressed as a two-parameter ansdtz for the function 6c'- 

6c = ra.m{6o,6i-z}; (28) 

it will be referred to as ansdtzl. The parameter 6i describes how ZEL underestimates (5i<3) or overestimates 
(5i>3) the collapse rate of a mass element. In the limiting case 6i=6o, the plain ZEL curve is just shifted to 
cross the 6c axis at 6o. The PS case is also recovered in the limit 6i — > co. Inserting Eq. (28) into Eq. (27), 
X((7) can be written as: 

I[a)= [ [ dX dY V{X, Y) + ^ exp (- ~/° ) / / dX dYV{X,Y)exp (—Z - -zA . (29) 
J Jz<A/<T *o V 2(7^ J J Jz>A/<T 2 y 
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where Z = z/a, and A = 6i — 6o- Fig. 2 shows some X((7) curves for 5i-values 2, 2.5, 3 and 3.5 (dashed lines). 

It is easy to obtain the limiting behaviours of Eq. (29). In the limit of large mass variances, which corresponds 
to low masses. A/it is small, and then the first integral gives negligible contribution, while the second one turns 
out to be a constant, whose value is X((7) — > 0.33 6i/6o- In the limit of small variances, the second integral can 
be calculated in the variables X and Z: 

The sign of the constant (95q — 55^) determines whether T{a') continues growing at small a or has a maximum 
and then falls down. The former behaviour is present when 6i < •\/9/55o — 2.27; in this case the exponential 
cutoff of the dF/da term, Eq. (26), becomes exp[—56l/18a^], instead of the PS exp[— 5q/2(72] one. In the 
6i > •\/9/55o — 2.27 case, the second integral of Eq (29) becomes negligible at large a, while the first one 
approaches unity (for large A/it that is just the integral of a normalized probability), so that X((7) ~ 1. It is 
remarkable that, while the collapse of rare events is always nearly spherical, it does not guarantee the PS mass 
function to be recovered at large masses. 

In Fig. 2 the limiting behaviour, Eq. (30), is reported (dotted lines); it describes very well the relevant 
small-(7 features of the X((7) curves. In this way it is possible to estimate the position of the maximum of X((7) 
in the 6i > •\/9/55o case, which is 



(TMAX ^ - 95o2)/i8, (31) 

and its height. 



'max ^ 3.66 — (32) 



In this ansdtzl, being the 6c curve flat around z = Q, perturbations with a small initial shear collapse exactly 
as spherical ones. As a consequence, the MF reduces to the PS one {I{o') — > 1), unless the contribution from 
non-spherically collapsing perturbations, which is given by the second integral, overcomes the one from the 
quasi-spherical clumps. This is not realistic, as it would mean that the largest collapsed structures are not 
roughly spherical, like rich clusters, but rather elongated ones, like collapsed filaments. 

It is interesting to analyze what happens if the 6c curve is not exactly flat around z = 0, a realistic possibility. 
Let us consider a second ansdtz for the function 6c'- 

6c{z) = 6q-6'^z (33) 

where 62 is the slope of the function 6c at small z (Fig. 1, dot-dashed lines, 62=Q.\, 0.2, 0.3); PS is recovered 
for 6'2=Q. This ansdtz will be referred to as ansdtz2. Fig. 3 shows the resulting X((7) curves for 6'2=Q.\, 0.2 and 
0.3 (dashed lines). Putting Eq. (33) into Eq. (27), and integrating in the variables X and Z, in the small a 
limit one obtains: 

60V5 6'2Hl ( 61 \ 

where ^ = A6'2 /{^ + 4^2'^). In this way the PS exponential cutoff becomes exp[— 6l{l - ^)/2(72], so that the 
effective value oi 60 is lowered by a factor \/l — ^. It is interesting to note that, even for vanishing but non-null 
52-values, the X((7) function diverges at small a instead of approaching 1. This could seem in contradiction with 
the fact that the collapse of rare event is nearly spherical. However, the physical quantity is not X((7) but the 
MF, and an exponential divergence of X((7) at small ir-values corresponds to an effective change of 5o, which is 
vanishingly small when 62 approaches zero. Again, the fact that rare events collapse in a quasi-spherical way 
does not assure the PS to be recovered at large masses, but this time the added objects are quasi-spherical. 

To give more realistic and dynamically motivated expressions for the function 6c{x,y), two different models 
are proposed. First, in the evolution equation for the density, Eq. (5), the non-local shear-term can be expressed 
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locally by means of ZEL (Eq. 11). Then Eq. (5) can be solved numerically to give an estimate of the collapse 
time; this procedure will be referred to as ZEL model. It is easy to see that, as long as ZEL predicts a slower 
collapse than the spherical case, ZEL shear has not a strong effect and then the density diverges nearly at 
the time predicted by the spherical model, while in the opposite case the divergence of ZEL shear drives the 
collapse, that will then happen nearly at the ZEL-predicted time. In practice, ZEL model is a dynamically 
motivated way to continuously join the spherical and ZEL regimes together. Inverting Eq. (19) it is possible to 
find the corresponding 6c function, which is shown in Fig. 1 (heavy dashed line); note however that this time 
6c is not simply a function oi z = x -\- 2y, except at large shears, but as a matter of fact the difference from 
the plain z-dependence is small; the curve in Fig. 1 is representative of the mean behaviour. As expected, ZEL 
model is very similar to the ansdtzl with 5i=3, with a non-vanishing slope at small x and y. The ZEL model 
X((7) curve is shown in Fig. 2 (heavy line); again it is very similar to the ansdtzl 5i=3 curve, apart from the 
weak exponential rise at small a, which takes place when the PS mass function has already been cut off. As a 
consequence, all the features of ansdtzl 5i=3 described before (Eqs. 31 and 32) can be applied to ZEL model 
to quantify its relevant features. 

The second model is the well-known homogeneous ellipsoid model (e.g., Peebles 1980); it will be called ELL 
model. Care has to be used in introducing it because, while a spherical perturbation can be consistently put in 
an otherwise homogeneous Universe, this is not possible for a triaxial overdensity. Again, what has to be done is 
to "extract" an ellipsoid from a whole perturbed Universe. Notably, the dynamical equations for the evolution 
of an ellipsoid are formally the same as the ones for the evolution of a general Lagrangian mass element (Peebles 
1980); what is totally different is the role of the potential, which is what gives informations about the rest of 
the world. Following Bond &: Myers (1993a), the gravitational potential can be expanded in a Taylor series: 

^(x) = ^0 + -^a;,- + ^ ^ f* XiXj + ... (35) 
axi 2 axidxj 

The first term is a constant and has no importance; the second term induces just a bulk motion on the mass 
element. The third term, which is a quadratic form and contains the Laplacian of the potential (its trace) and 
the tidal field (its traceless part), is the first and only relevant term which influences the internal properties 
of a mass element (no further derivatives of the potential enter the Lagrangian evolution equations of a mass 
element, Eqs. 1, 2 and 3). This is splitted in two parts: 



<l> = <I>EXT + <I>INT- (36) 

The first term, traceless, describes the external tides and is dominant in the early evolution; the other term, 
which dominates in the non-linear phase, describes internal dynamics, and can be written as the potential of 
a homogeneous ellipsoid, which is possible as it is a quadratic form. Initial conditions are set up so as to be 
consistent with ZEL at early times; again they are given by the three eigenvalues A,-. 

The equations of motion of the ellipsoid, which are reported in Appendix A, have been numerically integrated 
to determine the collapse time. Consistently with what has been done before, the collapse time has been defined 
as the instant at which the density becomes very large, which happens when the ellipsoid collapses along its 
shortest axis; this point will be discussed at the end of §4. The resulting 6c function is shown in Fig. 1 (heavy 
continuous line); again, the curve in Fig. 1 is just representative, as 6c is not simply a function of a; + 2y, but 
the difference from this behaviour is not strong. This time the slope of 6c at small x and y is larger than the 
previous cases; moreover. Ell model predicts slower large-shear collapses than the pure ZEL prediction, which 
is probably due to the strong non-radial motions which are important in these cases. Fig. 3 shows the resulting 
X((7) function (heavy line), in which the exponential rise at small variances is apparent; ansdtz2 curves with a 62 
parameter around 0.2 give a qualitatively good fit. As the asymptotical behaviour of ELL model at large x and 
y has turned out difficult to determine, and given that only the small-shear part is relevant for our purposes, 
plain ZEL prediction has been forced after x + 2y = i; this causes the small bump in the ELL X((7) curve at 
(7 ~ 4 and determines the asymptotical large-variance behaviour, which is the same as the ZEL model one. The 
pure exponential behaviour of ansdtz2 X((7) takes place at very small ir-values. It is then useful to compare the 
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ELL X((7) function directly to the one necessary to make the MF be equal to the PS one with an effective 60 
value , say 6e, instead of the 1.69 one: 



^W = ^expfM_ii). (37) 



60 V 20-2 

Fig. 3 shows the X((7) curves corresponding to 5b = 1.6, 1.5 and 1.4, values which have been found to fit 
several N-body results. The fit with the ELL curve is only rough, but considering that the small-ir difference can 
be lost in N-body simulations because of limited statistics, and that the it ~ 1 difference, which is slightly more 
than a factor 2, can be lowered by a correct normalization procedure (non-spherical objects are overcounted) 
and by fragmentation of clumps, it is reasonable to think that such an ELL MF could be roughly fit by a 5b < 60 
PS curve. 



4 Results and Discussion 

The last quantity needed to calculate the MF (Eq. 23) is the cosmological term da/dM; it introduces the 
dependence of the MF on the power spectrum of perturbations. It is useful to consider first scale-invariant 
power-law power spectra: P{k) oc k" . In this case: 

<7(M) = (M/Mo)-(^+")/^ (38) 
where Mq is the mass scale at which (7=1. Then our new MF takes the form: 

-U^l) (- i^f) ^ i^f) ^ <-) 

where = Mo(2/52)3/(3+n) expected, the corrected MF is scale invariant as long as the spectrum is scale 
invariant. Spectra have been normalized by the relation (7(8/i~^Mpc)=l (on a top-hat sphere) with /i=0.5, and 
gaussian smoothing is assumed. Fig. 4 shows the results for the PS (dashed), 5b = 1.5 PS (dotted), ZEL (dot- 
dashed) and ELL (continuous) MFs, for spectral indexes n from —2 to 1. The enhancement of the large-mass 
part is apparent in all the cases; in particular, the ZEL MF predicts more structures than the PS one around 
Mh, but reduces to it at large masses, while the ELL MF is has its large-mass tail shifted toward large masses, 
and is well-fitted in that part by the 5b = 1.5 curve. 

Because of Eq. (38), mass scales are more diluted for small spectral indexes, so that both the extent of 
the bump in the ZEL case and the shift of the large-mass part in the ELL case are larger for n = —2. With 
the aid of the two ansdtz presented in the previous section, it is possible to quantify in a more precise way 
the two features. Using ansdtzl with 61 = Z for the ZEL X((7) function, the maximum of the bump is at 
Mm AX = (5^/5)3/(3+") Mh. ~ 0.573/(3+") Mh. ^jjjj maximum height is a factor ~ 3.03 (Eq. 32). For 

the ELL mass function, it suffices to say that the shift of M* is (5o/5b)^^^^^"' (see Eq. 37), which is larger for 
smaller spectral indexes. 

All these features are recovered in more realistic power spectra. As an example, the standard cold dark 
matter (CDM) model has been used. The transfer function has been taken from Bardeen et al. (1986); f2 = 1, 
h = 0.5 and 6=1 have been assumed; gaussian smoothing has been used; the normalization has been fixed 
by aTH{^h~^Mpc) = 1, which is roughly consistent with COBE, given 6=1. Fig. 5 shows the four MFs 
(PS, 5b = 1.5 PS, ZEL and ELL) at different times; as expected, at earlier times, when the logarithmic slope 
of the spectrum at the collapsing scale is lower, the ZEL bump is slightly more diluted and the ELL shift of 
Mh. is slightly more pronounced. In any case, a rather firm conclusion can be drawn: because of the fact that 
non-spherical collapses are faster than spherical ones, more large-mass objects are present with respect to the 
PS mass function. This is consistent with the findings of Lucchin &: Matarrese (1988), who have found an 
increase of M* due to the non-Gaussianity of the evolved matter field, and with Cavaliere &: Menci (1994), who 
have found a similar effect by taking into account the coagulation of clumps. 
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The validity of the PS-like approach we have used is limited to the large-mass part of the MF (M > M^,) for 
a number of reasons. Single-stream dynamics has been assumed, which is valid if a mass element has a small 
probability of residing in multi-stream regions; Kofman et al. (1994) have shown that multi-stream regions 
are rare for mass variances smaller than it ~ 1, which means that a PS-like approach can be meaningful only 
for masses larger than . Furthermore, coagulation and fragmentation of clumps have been neglected. It is 
reasonable to assume that coagulation is negligible in single-stream regime (clumps have to come in touch to 
coagulate), while fragmentation of clumps is expected to be negligible when a mass element collapses in a quasi- 
spherical way, so that tides can not be so efficient to disrupt it. Another assumption is that the collapse time 
could be predicted from local initial conditions. Notably, a mass element spends most of its time before collapse 
in the linear and quasi-linear regime, where single-stream regions dominate and ZEL description is accurate, 
while the following highly non-linear phases are much faster. As a consequence, the collapse time predictions 
that have been presented here are likely to be realistic, while subsequent fragmentation and coagulation histories 
can give severe problems at small masses. A last assumption, already discussed before, is that PS normalization 
can be applied. This is surely wrong at small masses by a factor 1.84, but is probably accurate in the large-mass 
part, where the collapses are quasi-spherical. The application of excursion set theory can solve this last problem, 
but this time the diffusion of trajectories in the Ai, A2 and A3 (or 6, x and y) space has to be considered, and 
the whole 6c{x,y) function has to be used as an absorbing barrier. 

There is a number of points which deserve further discussion. On the actuality of Bertschinger &: Jain's 
collapse theorem, which states that the spherical perturbation is the slowest to collapse, there is not full 
agreement in literature. In practice, that theorem is valid for a fluid element, which has no structure by 
deflnition, and is then appropriate for our approach, in which the power spectrum has to be truncated to avoid 
strong non-linearities. Different reasonings, which can end up in different conclusions, apply to an extended 
structure. For instance, Peebles (1990) and Bond &: Myers (1993a) claim the shear to slow-down the collapse. 
The former's result will be discussed below; the latters, applying the ELL model to the collapse of density 
peaks to determine the instant at which full virialization is reached, have chosen to stop the evolution of every 
collapsing axis of the ellipsoid at a certain point, mimicking in this way a process of partial virialization, and 
to wait for the collapse of all the axes to deflne the collapse time. As a consequence, a non-spherical collapse 
is slower than the corresponding zero-shear spherical one. Actually, after flrst-axis collapse, because of partial 
virialization, matter can not be assumed any more cold and in single-stream regime, and then ELL model ought 
not to be extrapolated beyond that moment; nonetheless, predictions based on that model have been found in 
satisfactory agreement with N-body simulations (Bond &: Myers 1993b). 

In the present work, flrst-axis collapse has been chosen here as a suitable deflnition of collapse time. This 
different choice is mainly motivated by the fact that large overdensities and not fully virialized structures are 
searched. In practice, a mass element, as the ones that have been considered here, has no internal properties other 
than density, expansion and shear, so that the only meaningful deflnition of collapse is when those quantities 
reach large values; this happens at the moment corresponding to the flrst-axis collapse. It has been often said 
that such a collapse is only kinematical, in the sense that matter collides because of its initial velocity, while the 
internal potential does not diverge (Shandarin &: Zel'dovich 1989). On the other hand, the flrst-axis collapse is 
accompanied by shock waves in the baryonic component (Zel'dovich 1970), and after that moment matter can 
not be described any more as a cold fluid in single-stream regime, so that this is the true moment of transition 
to highly non-linear dynamics. Another objection which could be done is that in this way pancake-like objects 
are counted. It is to be stressed, however, that pancakes are only a transient feature of a self-gravitating system 
(Shandarin &: Zel'dovich 1989), which disappears after the collapse. So, a collapsing mass element will only 
show a fllamentary phase in the flrst instants of collapse, as is seen in N-body simulation of clusters (e.g.. White 
1994), but will then become a clumpy object; the subsequent possible fragmentation of the clump is the really 
worrying feature. 

A last assumption which has been made is that small-scale structure can be safely neglected by truncating the 
power spectrum at some wavelength. Taking small-scale structure into account, arguments based on least-action 
calculations (Peebles 1990) and dynamical friction (Antonuccio &: Colafrancesco 1993) lead to the conclusion 
that a collapsing structure is severely slowed down just after it has detached from Hubble expansion. Peebles 
(1990) has also shown that the presence of external tides induces strong non-radial flows in an otherwise spherical 
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overdensity, so that the collapse is slowed down and then non-spherical collapses are slower than spherical ones, 
a conclusion which is again opposite to Bertschinger &: Jain's collapse theorem. By the way, it is to be noted 
that Peebles' and Bertschinger &: Jain's statements have been obtained in quite different frameworks, so that a 
comparison of the two is not straightforward. Nonetheless, it is to be retained that, as more small-scale power 
is present, the collapse of a perturbation is slowed down in a way that could inhibit the positive effect of the 
shear. 

5 Summary and Conclusions 

In this paper the original Press &: Schechter (1974) approach to the mass function has been extended to 
non-spherical dynamics. Dynamical approximations and models of non-spherical collapse have been used to 
determine the instant of collapse of a Lagrangian mass element, starting from the local eigenvalues of the 
deformation tensor, which are the suitable initial conditions for local dynamical predictions. The Zel'dovich 
(ZEL) approximation, which is a true self-consistent first-order Lagrangian and local approximation, has been 
used as a guide for the determination of the collapse time. As ZEL does not reduce to the exact spherical model 
in the zero-shear case, two ansdtz for the true collapse time have been analysed in order to understand some 
relevant features. For a realistic estimate of the collapse time, two models have been used, one based on the 
numerical solution of the exact Lagrangian equation for the density contrast, but with the shear given by ZEL 
(ZEL model), the other based on the numerical solution of the homogeneous ellipsoid model up to the first-axis 
collapse (ELL model). Two fundamental merits of these ZEL inspired models, with respect to the spherical 
collapse, are the fact that collapsing clumps are consistently extracted from a perturbed field, instead of being 
put in an otherwise homogeneous Universe, and that the shear of the perturbation, which is the kinematical 
quantity that is missed by the spherical model, is taken into account. 

As a result, the collapse time of general non-spherical perturbations is shortened, in agreement with Bert- 
schinger &: Jain's collapse theorem, and then more large-mass objects are formed; ZEL model predicts more 
objects to be present around the knee of the MF, the ELL model, which is probably more realistic, is well 
fitted in the large-mass end by a PS function with a lower value of the effective 6c parameter, in the range 
1.4-^1.6. This last prediction is in agreement with many N-body simulations: e.g. Efstathiou &: Rees (1988) 
have found 5c = 1.33, Carlberg &; Couchman (1989) 1.44, Bond &; Myers (1993b) 1.58, Klypin et al. (1994) 1.33, 
while Efstathiou et al. (1988), Bond et al. (1991) and Lacey &: Cole (1994) have used the 1.69 spherical value. 
However, in all those cases the PS mass function has been treated as a phenomenological function, and no 
dynamical interpretation has never been given for 5c-values smaller than 1.69. 

When small-scale structure is present, dynamical friction and previrialization events can slow down the 
collapse rate of a perturbation, thus increasing the effective value of the 6c parameter. As a matter of fact, 
Klypin et al. (1994) have found their N-body MFs to be fitted at different times by PS functions with different 
5c-values, larger at later times. However, this result could be caused, at least in part, by their clump-finding 
algorithm, which can underestimate the number of large-mass objects. As mentioned in the introduction, also 
Jain &: Bertschinger (1993) have found their N-body MFs to give more large-mass objects than the pure PS 
prediction at large redshifts, consistently with their second-order calculation of the power spectrum, while at 
later times the PS curve has been found to give a good fit and even overestimate the N-body curve. In their 
case, the earlier-time N-body MFs are not well fitted by PS with a smaller 6c- It could be speculated then that 
the restoration of PS shape at later time is a consequence of the fact that previrialization tends to contrast the 
positive effect of the shear, so that at a certain moment spherical and non-spherical perturbations can collapse 
at the same time. 

To conclude, the PS mass function with 5c-values in the range 1.4-^1.6 is confirmed to give a good first 
approximation of the MF of large cosmological structures; besides, the use of non-spherical collapse models 
gives a dynamical explanation of the lower 6c parameters used to fit some N-body simulations, though the true 
shape of the mass function can be somehow different from the PS one. The possible decrease with time of the 
number of large-mass objects, or equivalently the increase of the effective 6c at later times, can have interesting 
consequences on the abundance of objects at high redshifts, as quasars (Efstathiou &: Rees 1988) or rich clusters 
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of galaxies (Evrard 1989; Peebles, Daly &: Juszkiewicz 1989). In particular, it could explain, at least in part, 
why rich clusters can be found at very high redshifts, which seems incompatible with their abundance in the 
near Universe. 



The author thanks Sabino Matarrese for his help and encouragement, and Stefano Borgani, Marco Bruni, 
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Appendix A: Equations for the Homogeneous Ellipsoid Model 

In this appendix the evolution equations for an homogeneous ellipsoid, ELL model, are given (for a derivation 
see Bond &: Myers 1993a). Let ri[t) = a,(t)5,- the physical i-th. coordinate of the outer surface of the ellipsoid, 
in the principal-axes system; q,- is the initial — Lagrangian — position and a,- is the expansion factor for the 
i-th axis. The universal scale factor is called a. Then the evolution equations for the a,- factors are: 

(40) 
where 

6 = - 1 (41) 

aia2a3 

b'i = ^[aiajakRoiaha'^j, 4)-!] i + j + h (42) 

A;. = -^ (^ + aoA.) (43) 

Here the A,- parameters are the eigenvalues of the initial deformation tensor, dui{ao) / dqi = diag (Ai, A2, A3), 
and Ro is the Carlson's elliptic integral: 



d ai 



-iwGgai 



1 ^ 

- H \- —6- 

3 3 2 



A' 
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(44) 



this integral has been calculated by means of the routine given by Press &: Teukolsky (1990). 
The initial conditions have been set up in order to be consistent with ZEL at early times: 



a,(to) = ao(l + aoA,) 



(45) 



ddi , , if da \ , , , 9 , ^ 



(46) 



Eq. (40) reduces to the usual Friedmann equation if all A,- = (which implies of course 5q = 0). The second 
term between squared parentheses in the rhs of Eq. (40) gives the spherical contribution from the overdensity, 
the third and fourth ones give the contribution from internal and external tides. This last term is evolved 
according to linear theory (Eq. 43), which corresponds to holding the peculiar potential constant; this is a good 
approximation as long as the perturbation is in the quasi-linear regime (e.g., Bagla &: Padmanabhan 1994), 
while in the highly non-linear regime the second and third term, the internal non-linear ones, dominate. Bond 
and Myers (1993a) have verified that changing this assumption with other reasonable ones does not affect the 
result. 

Figure captions 

Figure 1: 6c curves as a function oi z = x + 2y: plain ZEL prediction (dashed), ansdtzl curves (dotted), 
ansdtz2 curves (dot-dashed), ZEL model (heavy dashed) and ELL model (heavy continuous). 

Figure 2: X((7) curves of ansdtzl, with 6i=2, 2.5, 3 and 3.5 (dashed lines); dotted lines represent the 
corresponding limiting behaviour, Eq. (30). Heavy continuous line: ZEL model. 

Figure 3: X((7) curves of ansdtz2, with 6'2=Q.l, 0.2, and 0.3 (dashed lines). Continuous lines correspond to 
PS curves (Eq. 37) with 5b = 1.6 (left), 1.5 (center) and 1.4 (right). Heavy continuous line: ELL model. 

Figure 4: The mass function for power-law spectra [n = —2 and 0): PS (dashed), PS with 5b = 1.5 (dotted), 
ZEL model (dot-dashed) and ELL model (continuous). 

Figure 5: The mass function for a standard CDM spectrum, at different times: PS (dashed), PS with 
5b = 1.5 (dotted), ZEL model (dot-dashed) and ELL model (continuous). 
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